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ABSTRACT 



The problems of diffraction of an electromagnetic plane wave by a perfectly conducting 
circular disk and its complementary problem, diffraction by a circular hole in an infinite 
conducting plate, are rigorously solved using the method of the kobayashi petential. Th 
mathematical formulation involves dual integral equations derived from the potential in- 
tegral and boundary condition on the plane where a disk or hole is located. The field is 
expressed by a linear combination of function which satisfy the required boundary condi- 
tions except on the disk or hole. It may also be varified that the solution for the disk and 
the hole satisfy Babinet's principal. The weighting function in the potential integrals are 
determined by applying the properties of the Weber-Schafheitlin's integral and the solua- 
tion is obtained in the form of a matrix equation. Matrix elements of the equations for the 
expansion coefficients are given by three kinds of infinite integrals and the series soluation 
for these infinite integral are derived. For the varification of these series soluation, the 
numerical integral are also derived and the results are computed numerically for confor- 
mation, which is fairly good. Illustrative computations are given for the far diffracted field 
pattern. The results of the far-field patterns are compared with the results obtained from 
physical optics (PO). The agreement is fairly good. 
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Chapter 1 
Introduction 



Over several decades, electromagnetic (EM) scattering from a circular disk has attracted 
researchers with a strong intrest in the (monostatic or backscattering) radar cross section 
(RCS's) of large dynamic range. Usually, as reviewed by Duan and Rahmat-Samii [1], the 
following groups of techniques are applied to the analysis of the electromagnetic scattering 
from circular disk. The first type is the physical optics (PO) [2] which is an approximate 
technique and is accurate for predicting the far-filed pattern near the main beam. The 
second type is the physical theory of diffraction (PTD) [3] which is more accurate than the 
PO technique since the equivalent edge current is applied and the caustic singularities in the 
original ray tracing are eliminated. This method is further modified (a) by Ando [4] using 
equivalent edge currents, (b) by Mitzner [5] utilizing the incremental length diffraction 
coefficients, and (c) by Michaeli [6,7] using surface-to-edge integral and the fringe current 
radiation integral over the ray coordinates instead of the normal coordinates. The third 
type is the geometrical theory of diffraction (GTD) [9-13], which has similar accuracy to 
the PTD. This method was also modified into (a) uniform geometrical theory of diffraction 
by Kouyoumijian and Pathak [14], (b) uniform asymptotic technique by Ahluwalia et al. 
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[15] and Lee and Deschamps [16], and (c) High-Order Geometrical Theory of Diffraction by 
Bechtel [17] and Ryan and Peters [18] (of the 1st order), by Knott et al. [19] (of the second 
order), and by Marsland et al. [20] (of the higher-order). The fourth type is the method of 
moments (MoM) or moment method (MM) [1] that is considered to be numerically exact. 
The Hybrid Asymptotic Moment Method is implemented by Kim and Thiele first who 
found the induced currents on the scatterer surface. This method was further modified 
by Kaye, Murthy and Thiele [21,22], and thus the fifth type of methods was formed i.e., 
the hybrid-iterative method which employs the magnetic field integral equation for the 
induced currents to solve the scattering problem. The diffraction of electromagnetic 

waves by a circular disk and a circular hole is solved by Bouwkamp [23] by using the power 
series expansion, by Meixner and Anderjewski [24], by Andrejewski [25] and by Flammer 
[26] by using the spheroidal wave function, and by Levine and Schwinger [27] by using 
the variation method. The variation method is not sufficient since it depends on the trial 
functions. The spheroidal wave function is not adequate in extending it to a complicated 
problem, but give great numerical results for the case of circular disk. Following the idea 
of Meixner and Anderjewski, Nomura and Katsura [28] reformulated this problem by using 
the Weber-Schafheitlin's discontinuous integral (Kobayashi Potential). Here we discuss 
an alternative method to that by Nomura and Katsura. Most of the researcher have 

used an integral equation for unknown equivalent surface current density on the aperture 
or disk. This integral equation is reduce to matrix equation via the method of moments 
(MoM). In this dissertation, rigorous solution to the problem of a plane wave scattering by 
a circular conducting disk and its complementary problem, diffraction by a circular hole 
on a perfectly conducting plane, are derived using the method of the Kobayashi potential 
method [29], [30]. This method has been applied to various kinds of problem such as the 
potential problems of electrified circular disks [31], [32], the diffraction of acoustic waves by 
a circular disk [33], the diffraction of electromagnetic plane wave by a rectangular plate and 
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a rectangular hole in conducting plate [34] , and the diffraction of acoustic plane wave by a 
rectangular plate [35], [36]. The Kobayashi potential has also been used for the diffraction 
of electromagnetic waves by a thick slit [37], a flanged parallel-plate waveguide [38], an N- 
slit array [39], etc. The Kobayashi potential method resembles the MoM in its spectrum 

domain, but the formulation is different. The MoM is based on an integral equation, 
whereas the Kobayashi potential method starts from the dual integral equation. The MoM 
in a space domain has been used mostly in the diffraction problems of electromagnetic 
waves. We can cite the following advantages of the Kobayashi potential method over the 
current numerical techniques (mainly over MoM). 

1. In contrast to the MoM in a space domain, the Kobayashi potential method does not 
involve singularities of the Green's functions, so we can obtain very accurate results. 

2. Since each function involved in the integrand of the potential functions satisfies a part 
of the required boundary condition, the convergence is very rapid. In this respect, 
the present method may be regarded as eigenfunction expansion of the geometries. 

3. As in two-dimensional (2-D) problem, the Kobayashi potential methods may be ap- 
plied to more complicated problem with related configurations. These problems may 
be formulated in a manner similar to the eigenfunction expansions in cylindrical and 
spherical geometries. 

4. For 2-D problems, the solution to a two-slit diffraction can be used to predict the cou- 
pling between the slits asymptotically [40] . This is also expected in three-dimensional 
(3-D) problems. 

The disadvantage is that, the tractable geometries of this method are limited to special 
shapes like rectangular and circular plates and their related geometries. A similar situation 
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is seen for other conventional eigenfunction expansions. 

In this dissertation, in which we discuss the diffraction of electromagnetic waves by a 
circular disk and circular hole, the solution begins by introducing the Fourier sine and cosine 
transforms of the tangential components of the vector potentials. From the requirement 
of the boundary conditions on the plane exterior to the disk or hole, we obtained the dual 
integral equations for the transformed functions (or weighting function). The equation are 
solved by using the properties of the Weber-Schafheitlin discontinuous integrals. At this 
step, we can incorporate the required edge condition into the solution. The results include 
two kinds of arbitrary discrete parameters, so that the general solution is obtained by 
superposing these results. By imposing the remaining boundary conditions on the circular 
disk or on the circular hole, we have a matrix equation for the expansion coefficients. 
Matrix elements are given by an infinite integral, which are then expanded into infinite 
series, which is more convenient for numerical computation. The numerical results has 
been presented for the far-filed pattern and compared the results of the far-filed pattern 
with the corresponding physical optics (PO) method. More detailed and the mathematical 
formulation will be examined in the next chapter. 
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Chapter 2 

Diffraction of Electromagnetic Wave 
by Circular Disk and Circular Hole 

The diffraction of electromagnetic plane wave by a circular disk and its complementary 
problem, diffraction by circular hole in a perfectly conducting plane, are classical problems 
and has been formulated rigorously by using eigenfunction (spheroidal function) expansion 
[24]~[26] and the method of the Kobayashi potential [28]. Meixner showed how to formulate 
the problem as the boundary value problem by using the rectangular components of Hertz 
vector. He first split the Hertz vector into two parts. One is associated with the incident 
wave and other is related to the scattered wave. He first noticed that the field on the circular 
aperture can be expressed in term of two dimensional scalar wave function and expanded 
this function in term of the spheroidal function and he expressed it by taking into account 
of the assumed expression on the aperture. The auxiliary function representing aperture 
field was determined by imposing the edge condition. By using this formula Meixner and 
Andrewski gave the numerical results. 
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Following this idea Nomura and Katsura [28] reformulated this problem by using the 
Weber-Schafheitlin's discontinuous integrals (Kobayashi Potential). We discuss here an 
alternative method to that by Nomura and Katsura. We have the longitudinal components 
of the vector potentials of the electric and magnetic type and drives the dual integral 
equations for the surface field on the plane where disk or aperture is located. The solution 
for one of the pair equations, is expanded in terms of the functions which is derived by taking 
into account the Maxwell's equations, discontinuous properties of the Weber-Schafheitlin's 
integral and the required edge conditions. The corresponding spectral functions of the 
surface field may be derived by applying the vector Hankel transform discussed by Chew 
and Kong [41], [42]. This determines the weighting functions of the vector potentials. The 
expansion coefficients are determined from the another of the dual equation by using the 
projection of the function space. Then we can reduce the problem into matrix equation. 
Matrix elements are given by infinite integrals and can be expressed in term of infinite series 
which is convenient for numerical computation. Since the formulation is more compact than 
the method by Nomura and Katsura, it is promising to apply to more complex problem. 



2.1 Statement of the problem 



The geometry of the problem and the associated coordinates are described in Fig.2.1, 
where radius of the hole and disk is a and thickness of the conducting plane is assumed 
to be negligibly small. There are two kind of polarizations for the incident plane wave. 
Since we consider both scattering problem of disk and complementary circular aperture, 
we list the electromagnetic filed components and the corresponding vector potentials A z 
and F z for the incident wave and wave reflected from the infinite conducting plane located 
at z — 0. The relation between (E,H) and (A z , F z ) are given by 
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H = -VxA (2.1) 
E = --V x F (2.2) 

Since we are using the longitudinal components of the vector potentials of the electric and 
magnetics type, hence the above equations can be rewritten as 

l,dA z . dA z , l f dA z . dA z , 
l f dF z . dF z .^ l f dF z . dF z .^ 

E = ~(V as "^ = ~W""V^ ( } 

These relation can be used to find the corresponding electromagnetic field from the vector 
potentials. 



2.2 Plane Waves in Term of Cylindrical Wave Func- 
tions 

To facilitate to enforce the boundary condition we express the plane wave in the circular 
coordinates by using the generating function for Bessel function, given by 



exp 



1 oo 

■(*-t)]= E Ux)t m (2.5) 

~" m=— oo 

Putting t = exp(j| — j<f>) in above equation, we can write 

oo 

exp[jxcos0] = J2 j m J m (x)exp[-jm(f)} (2.6) 

m=— oo 

Using the recurrence relation for the Bessel function, that is 

J^ m {x) = (-l) m J m (x) 
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Eq.(2.6) can be written as 

oo 

exp[j/tpsin#o c °s0] = X! j m Jm(i^p sin 6> ) exp(—jm<f>) 

m=— oo 

oo 

= ^2 e m j m J m (K,p sin 6 )cosm<i) (2.7) 

m=0 

where e m is Neumann's constant defined by e m = 1 for m = and e m = 2 for m > 1, and 
the symbol # denotes the polar angle of the incident wave. 



2.3 Perpendicular (Horizontal or E) Polarization 

Let us assume that the plane of incidence lies in xz plane, then the incident electromag- 
netic plane wave on perfectly conducting plane are given by 

E* = i y Ei exp[j kx sin 9o + j kz cos 9 ] (2.8) 

EE = Y Ei(i x cos9 — i z sin 6 )exp[jK,x sin 9 + jn,z cos 6 ] (2.9) 

Similarly the reflected fields can be expressed as 

E r = — % y E\ exp[jnx sin 9 — jnz cos 9o] (2-10) 

EE = Y Ei(i x cos 9 + i z sin 9 ) exp[jnx sin 9 — jnz cos 9 ] (2-11) 

Using Eq.(2.3) and (2.4), the corresponding vector potential can written as 

€qE\ 

Ft = — ex^Unx sin 9q+ jnz cos 9q] (2.12) 

jKsm9o 

€qE\ 

F r = — expljnx sin 9q — jnz cos 9 ] (2.13) 

]Ksm9 

Using the transformation from rectangular to cylindrical coordinates, the tangential com- 
ponents of the electromagnetic field on the plane z = are given by 

E p = E y sin<f> = ^isin^exp [jnp cos <p sin 9 ] (2-14) 



H p = H x cos = Y E 1 cos #o cos exp [jnp cos sin 9 ] 
= E y cos 4> = Ei cos exp [j up cos sin 9 ] 
= —H x sin = — Yo-Ei cos #0 cos exp[jKp cos sin O ] 

From the above equation, it can readily carried out 

H p = YqE^ cos 9 
= -Y E P cos 9 

Now applying the wave transformation, Eq.(2.7), we have 



(2.15) 
(2.16) 
(2.17) 

(2.18) 
(2.19) 



m 



h; = y e; cos 9 

00 

= -jY Ei cos6» J2 £ m j m J' m ( K P sin0 o)cos 

m=0 

H r 6 = -Y El cos 9 



1716 



m 



= jY Eicos9 V e m j m — ^— r J m (/tp sin 9 ) sin mc, 
— np sin 6*0 



m=0 



where we have used the relation, 



X / 1\ rX / 1\ 

-(l + -)exp[-( t --) 

1/ 1\ rX / 1\ 

5 -; exp k '-7 



= 51 rnJ m (x)t r ' 

m=— 00 

00 

= E ■4(*)* m 



(2.20) 
(2.21) 

(2.22) 
(2.23) 



where J' m {x) is the derivative with respect to the argument and Y is the intrinsic admittance 
of free space. 



2.4 Parallel (Vertical or H) Polarization 



In this case, the corresponding incident and reflected waves are given by 
E* = E 2 (i x cos9o — i z sin 9o)exp\jKX sin 9q+ jnz cos 9o\ 
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(2.24) 



H* = — Y E 2 i y exp\j kx sin 9 + jK,zcos9 ] (2.25) 

E r = E 2 (—i x cos #o — iz si n #o) exp[jnx sin 9 — jkz cos 9 ] (2.26) 

H r = — Y E 2 i y exp[jnx sin 9 — jnz cos 9o] (2.27) 

= ^°^ ^? expfmx sin fl + jkz cos fl ] (2.28) 

^ = ^°^ ^ 2 exp[j/ta;sin^o — j/t^cos^o] (2.29) 

J /t Sill C7g 

and the cylindrical coordinate expression of these waves are 

00 m 

H; = H;=jY E 2 J2^mJ m — J m («:psin^ o )sinm0 (2.30) 

m=0 Kp SmV 

E; = Z cos 9 H; (2.31) 

oo 

HI = H^ = jY Q E 2 ^€^ m J^Kpsme )oo8rn<f> (2.32) 

E; = -Z cos 9 H'; (2.33) 



2.5 The Expression for the Fields Scattered by a Disk 

We now discuss about our analytical method for predicting the filed scattered by a 
perfectly conducting disk on the plane at z — 0. We studied this problem by applying the 
Kobayashi potential [45]. 
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(a) 



(b) 



Figure 2.1: Diffraction of Plane Wave by a Circular Disk and a Circular Hole 

2.6 Spectrum Function of the Current Density on the 
Disk 

We assume the vector potential corresponding to the diffracted field as the superposition 
of the elementary solution for the wave function in the form 



^(p,<M) = ±A*oare*o / fcm(0 cosm0 + f sm (0 sinm0 
xJ m (PaO exp[=F^/f 2 -K 2 z ]f -1 df 

F?(p, 0, z) = e aJ2 9cm(0 c °s m0 + <7«n(f ) sin m. 

x J m (poO exp[=F^ 2 - K 2 2 a ]f _1 df 



(2.34) 



(2.35) 



where the upper and lower signs refer to the region z > and 2 < 0, respectively, and 
p a = ^ and £ a = ^ are the normalized variables with respect to the radius a of the disk. 
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Also /(£) and g(£) are the unknown spectrum functions and they are to be determined so 
that they satisfy all the required boundary condition. Eq.(2.34) and (2.35) are of the form 
of the Hankel transform for z — 0. It is seen that the tangential components of the electric 
filed derive from Eq.(2.34) and (2.35) satisfy the continuity on the plane z = due to the 
signs in front of Eq.(2.34). First we consider the surface field at the plane z = to derive 
the dual integral equations associated with them. By using the relation between the vector 
potentials and the electromagnetic field, give by 

(2.36) 



E = -juA-j— V(V- A) - -V x F 

ojjie e 

H = *V x A-juF- j — V(V-F) 
fx cope 



(2.37) 



E d A P , 0, 0) 



Substituting Eq.(2.34) and (2.35) into the above equation, the electromagnetic field com- 
ponents tangential to z— plane becomes 

oo 

\ E pc,m{Pa) cosm0 + E ps , m (p a ) sin 777 

m=0 

°° roc I r ^ „ i 

= jJ2 V^ 2 ~~ k2 [fcm(0 cosm0 + / sm (£) sinm^J J' m {p a O d i 

°° f°° r 1 777 

~H -9cm(xi)smm(f) + g sm (^)cosm(f) —J m (p a ^)d^ (2.38) 

oo 

{Etpcmipa) COS 7770 + £^ >m (p a ) SU1 777 

m=0 

°° poo , — 



t2 _ K 2 



i=0 

°° roc r 

+ J2 9cm(0 cos m(f) + g sm sin 7770 J' m (p a O d Z. 

oo 

H t(p, 0> °) = Y [ H pc,m(pa) cos 777 + H pSyTn sin 777 



"/cm(0 Sin 7710 + fsm COS 7770 — J m (paO d £ 



ipa 



(2.39) 



m=0 



= ±«Vb H / -fcm(Q Sin 7770 + / sm COS 7770 — J m (p a O d C 



^Pa 



K m=0 JO 
Xj'm(paO d t 



£ 2 - 0cro(f ) COS 7770 + g sm {0 SHI 777^ 



(2.40) 
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H$(p,(/),0) = [Hfamipa) cos m<f> + Hfarnipa) sin m(/) 



m=0 



m=0 



TkY / /cm(0 c °s m <P + Ism sin m<j> J' m (pai)di 



K m =0 J0 
£Pa 



'9cm{C) sinm0 + g sm (0 cosm0 



(2.41) 



In the above equations the upper and lower signs denote the values at z = + and z = 0_ , 
respectively. The tangential electric filed components take the same values at z = + and 
z =_. We rewrite Eq.(2.38)~(2.41) by using the matrices, given by 



Ep C ,m(Pa) 

E(f>s,m (Pa) 
-^ps,m(Pa) 

E(/)c,m{Pa) 



- f 

JO 



K p C) m (Pa) 
K<p s , m {pa) 



K ps,m(Pa) 
K<j)C,m(Pa) 



f C 

= 2Y 

Jo 



= 2Yr 



»/ 

Jo 

roo 

Jo 



H-{iPa) 
H + {i Pa ) 

~H-{ip a ) 

H-(tPa) 

H + {iPa) 
H + {iPa 



(2.42) 
(2.43) 



jVF^UO(^) -1 



(2.44) 



(2.45) 



where we use the current density K instead of surface magnetic field by using the relation 

are defined by 



K p = —1H^ and K$ = 2H p . The kernel matrix |^/ + (£p a )J and H (£p a ) 

J'j£ Pa ) ±§- a Jm{tPa) 

The required boundary conditions state that the current densities on the plane z = 
are zero for p a > 1 and the tangential components of the total electric field vanish on the 



(2.46) 
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disk. These are written as 



[CO r 

/ [h-(M 

JO 1 

H + {iPa) 



ps,m (0 
K<t>c,m(0 



Pa > 1 
Pa > 1 



(2.47) 
(2.48) 



(Pa) 



pc,m 



Zo 

2 j. 



H-^Pa) 



ps,m 



(Pa) 
(Pa) 



+ 

2 
+ 



E l pc , m (Pa) 
E <t>s,m(Paj 



roc r 
/ \H + (tPa) 

Jo L 



E ps ™(Pa) 



E\ 







Pa < i 



(2.49) 







Pa < 1 



(2.50) 



where E l pcm and -E* Sim denote the cosm0 and sinm0 parts of the incident wave E l , respec- 
tively, and the same is true for E l ^ c m and E l ^ s ^ m . The expression for theses factor can be 
obtained, and is given by 



E lc(pa) = -jE 2 cos 9 e m j m J' m (Kp a sin O ) 



EL(p a ) = -jE.emf 



m 



H Kp a Sill 9 

KM = -j J Eie m j m J^(/tp a sin0 o ) 



J m (/tp a sin6> ) 



E ls(Pa) = 3E2 cos 9 e m f 



m 



-Jm(Kp a Sm O ) 



(2.51) 
(2.52) 
(2.53) 
(2.54) 



np a sin 9 

Eq.(2.47)~(2.50) are the dual integral equations to determine the spectrum functions. 
To solve these equations, we expand K(p a ) by the function which satisfy the Maxwell's 
equations and the edge conditions. These functions can be found by taking into account 
the discontinuity property of the Weber-Schafheitlin's integrals. Once the expression for 
K(p a ) are established, the corresponding spectrum functions can be derived by applying 
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the vector Hankel transform introduced by Chew and Kong [41], [42]. Using the vector 
Hankel transform, [see Appendix A], Eq.(2.44) and (2.45) can be rewritten as 



K Pc , m (0 
^Q>s,m(0 



= f 

Jo 



Kpc t m (Pa) 
K^s^Pa) 
Kps t m (Pa) 

K(f>c,m(pa) 



Padpa 



Padpa 



(2.55) 



(2.56) 



It is noted that (K p , K$) satisfy the vector Helmholtz equation V 2 K + k 2 K. = in 
circular cylindrical coordinates, since K and H are relate by K = n x H on the plane 



0. Furthermore [K p ^K^) have the property K p ~ (1 — p^)a and ~ (1 



pI) 



near the edge of the disk. By taking into these facts, we substituting the Kernel matrix in 
Eq.(2.55), we can write 



J'JtPa) 



Jm{ipa) 



m 



. JUL 



Jm{£pa) 
JL(tPa) 



-^/9C,m(0 
^s,m(0 



S<% = 



(2.57) 



Using the recurrence relation for Bessel function, given by 



Tfl 

Jm-l(x) + J m+ l(x) = —J m (x) 

X 



and 



4-iW - J m +i(x) = J' m (x) 

and the properties of Weber Schafheitlin's discontinuous integral, [see Appendix B] , 
Eq.(2.47)and (2.48) can be satisfied if we choose, 

oo 

K pc , m (p a ) = T,[ A m nF~M ~ B mn Gi n (p a )] (2.58) 

n=0 

oo 

K ps , m (Pa) = Y,[ C ™nF mn (p a ) + B mn G+ m (p a )} (2.59) 

n=0 

oo 

K^ m (pa) = Y,[- A mnF+ n (p a )+B mn G mn (p a )} (2.60) 
n=0 

oo 

K^UPa) = Y,[ C rnnF+ n (p a ) + D mn G7 mn ( Pa )} (2.61) 



n=0 
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where 



,{Pa) = / 
JO 



J\m-l\(VPa)J\m-l\+2n+^(v) ± Jm+l(VPa) J m +2n+l(v)\v* df] (2.62) 
J\m-l\(VPa)J\ m -l\+2n+^(v) ± (VPa) ^m+2n+f fa)] »T**7 ( 2 - 63 ) 



It may readily be verified that F± n (p a ) = G± n (p a ) = for p a > 1, and F+ n (p a ) ~ 
(! - PaH, ~ (1 - Pa) ^ G+„(p a ) ~ (1 - § an d G mn {p a ) ~ (1 - near the 

edge p a ~ 1. Thus the above expressions satisfy one part of the dual integral equations 
with the unknown expansion coefficients A mn ~ D mn . To derive the spectrum functions 
/(£) and of the vector potentials, we first determine the spectrum functions of the 
current densities, since they are related to each other. Now substituting Eq.(2.58)~(2.61) 
into Eq.(2.55) and (2.56) and perform the integration, then the spectrum function of the 
current density is determined. The results are 



n=0 

oo 

n=0 

oo 

n=0 

oo 



n=0 



(2.64) 

(2.65) 
(2.66) 
(2.67) 



for m > 1 and 



k pc M = 2]>> „j 2n+| (Or i 

n=0 
oo 

K<f> s ,o(0 = 2]TA) m J 2n+| (£)r i 

n=0 

oo 

K P s,o(0 = -2£A,n</ 2 „ + |(£)r f 

n=0 

oo 

K^) = -2 £0,^+3(00 



(2.68) 
(2.69) 
(2.70) 
(2.71) 



n=0 
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for m = 0. In the above equations the functions S^j n (£) and r^ n (£) are defined by 

ZL(0 = [jm + 2n-\{0 ± J m +*+%[C k (2-72) 
r mn(0 = [j m +2n+\(Q ± J m+2n+§ (0] £~ 5 ( 2 -73) 

In deriving Eq.(2.64)~(2.67) and Eq.(2.68)~(2.71), we used the formula of the Hankel 
transform given by 

/ J m (ap)J m ((5p)pdp = — — — 
Jo a 

and the property of the delta function given by 



/ f(x)6(x-x')dx = f(x') 
Jo 



From Eq.(2.44) and (2.45), the spectral functions / cm (£) ~ 9sm(Q can be expressed in 
term of K pcs (£) ~ K^g). 

f cm (0 = ^K P cm(m (2-74) 
fsm(0 = ^K ps , m (m (2-75) 
fiUO = — =5=^^ m (0e (2.76) 
flUO = ^L^ /Wflg (2-77) 

These expression gives the relation between the weight function and the spectrum function 
of the current density which can be used for the derivation of the expansion coefficient. 



2.7 Derivation of the Expansion Coefficients 

The equation for the expansion coefficients can be obtained by applying the remaining 
boundary condition that the tangential components of the electric field vanish on the disk, 
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which is given by Eq.(2.49) and (2.50). By substituting Eq.(2.74)~(2.77) into Eq.(2.49) 
and (2.50) , we have the relation given by 

jy /£2_ K 2 ~ 
I »• n 



E Lm(Pa) 



^0 

2 j {) 



E ps,m(Pa) 



+ 

2 
+ 



Pa<l 



(2.78) 



K<Pc,m(0 



j\/f 2 -K 

= Pa < 1 



E ' P c, m {pa) 

E )>s,m{Pa) . 

roc r 

/ \H + (ZPa) 
JO L 

E ps,m(Pa) 

E \c,m\Pa) . 

Eq.(2.78) and (2.79 ) are projected into function space with element v™(p 2 a ) for E p and 
u™(pD for E^, [see Appendix c], then we obtain the matrix equations for the expansion 
coefficients A mn ~ D mn . The results are given by @finalout 



(2.79) 



I final out ^ A Tl 



n=0 



K JO 



J m+2n-\{£) + J m+2n+%{£) ^ J m+2p+\{0 ~ ( a p + 3 ) J m+2p+f (0 ^ 



T 7777, Z" 00 / r -I r -I 

T, A mn\ — / ~ ^ [ J m+2n-\ (0 + ^m+2n+§ (0 [J m +2 P -\ (0 ~ 4>+2p+§ (0 ^ 

n=0 ru JO - 

f°° 1 r i r i 1 

00 f jm f°° ^/£ 2 — /t 2 r 

+ 51 5 ™« i — / 7 J m+ 2 "+5 ^ ~ J ™+2n+f (0 

n=0 { JO C, 

-j* / OO ev ^ ? [^ + 2n + l(0 + ^ m+ 2n + |(0] [<^ + 2p-|(0 " (< + 1) J m+ 2 P+ I (0] « } 

= -j4Y E 2 cos 6 j m m J m+2p _i (/tsin6» ) + J m+2 p+|( Ksin6, o) (ftsinflo) - ^ 



J m+2p -l(0 + 4h-2 P +§(0J*; 



A) + J m+ 2p + |(/tsin6'o)](/tsin6'o) 2 
"i = l,2,3,---, p = 0, 1,2,3, 

(2.80) 



oo 

E 

n=0 



j_ r 

K JO 



J m+ 2n-h(0 + J , 



m+2n+S 



(0 <J m+ 2 P+ i (0 - « + 3) J m+2p+f (0 rf£ 
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poo 



1 r n 

00 f j roc JCl _ K 2 

+ E Ann I ^ X g2 [-Wan+lCO " ^ + 2n + |(0] [< " « + 3) J m+2p+| (£)] « 

/■oo 1 . IP , \ 

—jKmj^ C2^J C2 _ ^m+2n+|(0 + ^m+2n+|(0 ^m+2p+| (0 + ^m+2p+| (0 ^ J" 



n=0 



3 

= j4y ^ij m m[^ m+2 p+i(^sin6'o) + J m+ ^£(«sin0 o )J(Ksin0 o )~ 5 

m = l,2,3, p = 0,1,2,3, 

(2.81) 

f jffi roc I ,- , _ , 

E C mn \ — / ^ 2 ~ K 2 J m+2 n-| (0 + ^ m+ 2n+f (0 ^m+2p-i(0 + ^ m +2p+§ (0 ^ 
n=0 

/■oo 1 r , r , ~| 

J y/C2_ K 2 [ J m+2n-f(0 ~ ^ m+ 2n+|(0j [<^ m+ 2p-|(0 " « + 1) J m+2p+ | (0J ^| 

+ ^ A™ S — / -2 ^ m +2n+± (0 - ^ m +2n+f (0 J m+2p-i (0 + J m+2p+'l (0 ^ 

n=0 I > 



/"OO 1 r 

~i K J C ^J C2 _ ft f ^m+2n+±(£) + ^m+2n+§(0 



*™ J m+2p-§(0 - ("p" + 1 )^m+2p+|(0]^| 

= j4Yo£ij m [a^ m+2p -i(«sin0o) - « + 1) J m+2p+ |(/t sin# )] (/tsintfo)^ 

m = l,2,3, p = 0, 1,2,3, 



(2.82) 



1 /• 
E ^o™ - / 



1 ,00^ 



-J2n + §(0[-P-Vi(0 + (P+ l-5)J 2j9+ |(0l^ 



r "13 

= cos -pJ 2p+ i(/tsin6 l o) + (p + 1.5) J 2p+ s (/tsin^o) (/tsin# )~^ 



E ^0n K 



n=0 



o ^==^4 + i(0l-?Vi(0 + (p + i.5)j 2p+| (OJde 
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Y E 1 -pJ 2 p _i(k sin 9 Q ) + (p + 0.5) J 2 + s (resin 6 ) (resin6> ) 2 



P = 01,2,3,-.-, 



(2.83) 



where a™ = m + 2p. From Eq.(2.80)~(2.84) we find that all of the matrix elements have 
the form K(a, f3\) or G(a, /3, A), which are defined by 



These integral converge when a + (3 > — 1 and A > 1 for K(a, (3, A) and a + (3 > A — 1 
and A > — 1 for G(a, (3, A). By using the recurrence relation for the Bessel function, we can 
find that all the matrix elements contained in Eq.(2.80)~(2.84) converge for all indices. It 
is worthwhile to note that K(a, (3, A) and G(a, (3, A) can be integrated into infinite series 
which are more convenient for numerical computation. 

2.8 Far Field Expression 

The expression for the far filed can be drive by two different way. One method is to eval- 
uate the field radiated from the current density induced on the disk and the second method 
is to evaluate the expression of the vector potentials given in Eq.(2.34) and (2.35), directly 
by applying the stationary phase method of integration. Let the rectangular components 
of the current density be K x (p' a1 (p') and K y (p' a , <fi'), then the vector potential produce by 
this current density is given by 




(2.85) 



(2.84) 



A 




) 



(2.86) 



X 
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y = ~^k!s yi * PaA) r 



(2.87) 



where r is the radial distance between the observation and source points. 

Since we are interested in the far field, it can be shown that the radial distance r from 
any point on the source or scatterer to the observation point can be assumed to be parallel 
to the radial distance R from the origin to the observation point. In such cases the relation 
between the magnitude of r and R , given by 



R 2 - 2Rp' sin 9 cos(0 - 0' + p' 2 ) 
can be approximated, most commonly by [43] 

R — p' sin 9 cos(0 — <fi') for phase variations 
R for amplitude variation 

Substituting the above equation into Eq.(2.87)and (2.88), we have 

A x = p G (R) f K x (p' a ,<f)')exp[jKp' a sm9cos((j)-(j)')}dS' 
J s 

Ay = pqG (R) / K y (p' a ,(f)')exp[jKp' a sm9cos((f) - <p')]dS' 
J s 

where G (R) is defined by 



(2.88) 
(2.89) 



G Q (R) 



exp(—jK,R) 
AnR 



Using the transformation from cylindrical-to-rectangular coordinates, the rectangular com- 
ponents of the current density can be express as 



Since, 



K x{p' a A') = K p cos(f/ - K $ sin. <ft 
K v {PaA') = K p sin 0' + cos 

oo 

= [ K pc,m cos m(f)' + K pStJn sin m<j)' 

m=0 

oo 

[K<l>c,m cos m(f)' + Kfam sin m(j)' 
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m=0 



(2.90) 
(2.91) 

(2.92) 
(2.93) 



Hence, 

oo 

K x (p' a , 0') = [( K pc,m cos m<f>' + K pStrn sin m(f>) cos 0' - {K^ m cos ra0' 
+K ( p Sim sinm0') sin0' 

oo 

^(pI) 0') = 53 [(^pc,m cos m4> + K p . %m sin m0') sin 0' + (fQ, Cjm cos m0' 

m=0 

sinm0') cos0'] 

Substituting Eq.(2.91) and (2.92) into the above equation, we can write 

^(p, 0') = - ^ {(X>c,m - ifya.m) COs[(m - 1)0'] + (X pc , m + ify,, m ) 
2 m=0 

x cos[(m + 1)0'] 

+ ^ C ,m) sin[(m - 1)0'] + (x ps , m - fQ, c , m ) 
x sin[(m + 1)0']} 

J oo 

#y(p'» 0') = ^ E {-(^pc,m - ifya,m) Sm[(m - 1)0'] + (K pc ^ m + ^ s , m ) 
^ m=0 

x sin[(m + 1)0'] 

+ (K ps + ^Q> c ,m) cos[(m - 1)0'] - (if pSjm - i*Q, c , m ) 
x cos[(m + 1)0']} 

where we have use the relation, 

cos a cos = - cos(o; — /?) + - cos(a + (3) 

sin a sin /3 = ^ cos(o; — j3) — ^ cos(o; + /3) 

sin ct cos /3 = - sin(a; — /?) + - sin(o; + (3) 
Now substituting Eq.(2.97) and (2.98) into Eq.(2.89) and (2.90), we can write 

A x = ira 2 iioG (R) ]T j™" 1 / ^ 
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x {(#pcm ~ K <f>s,m)J m -i(Kp' a sin 0) cos[(m - 1)0] - (X pc , m + i^ s , m ) 
x J m+1 (np' a sin 6 1 ) cos[(m + 1)0] 

,171 + -fQ» c ,m) J m -i(/c/OoSin0) sin[(ra - 1)0] - {K ps ^ m - 
xj m+1 (/tp / a sin^)sin[(m + 1)0]} (2.98) 

Ay = TraVoCo^) E i m_1 / /W« 

m=0 J0 

x {-(K pCtm - J m _i(np' a sin 0) sin[(m - 1)0] - (K pc>m + i^ s , m ) 

x J m+1 (K,p' a sin 0) sin[(m + 1)0] 

,rn + ^Q. c ,m) J m -i(Kp' a sm0) cos[(m - 1)0] + (K pSyTn - l^ c , m ) 
x J m+1 (/tp' a sin0)cos[(m + 1)0]} (2.99) 

In deriving Eq.(2.99)and (2.100), we use the formula of the integral representation of the 
Bessel function, given by 



J n ( x ) — / exp(jx cos 9 + jn9)d9 

2ir Jo 



7T JO 



exp(jx cos 9) cos(n9)d9 



By applying the relation 



A e = A x cos 9 cos + Ay cos 9 sin 0, A^ = —A x sin + A y cos 
the spherical coordinate components of the vector potential becomes 

00 ra 

A e = 7ra 2 ^G (R)cos9j2j m - 1 pWa 

m=0 J ° 

x { (K pc . m - K (j>s ^J m - 1 (Kp' a sva9) - (K pc , m + K (j)S ^J m+1 (Kp' a sin9) 
x cos(m0) 

+ [(K ps , m + K^Jm-lWa^O) - (K ps 

,m K(j)C,rrij J-m+l {^Pa ^)] 

x sin(m0)} (2.100) 

00 /-a 

^ = na 2 n G (R)cos9j2j m - 1 pWa 

x {-[(-fr pc , m - i^ s , m ) J m _i(Kp' a sin0) + (K pc )j m+1 (Kp' a sin 9) 
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x sin(m0) 

+ [(K pS)m + if^, m ) J m _i(«Pa sin 6) + (K ps , m - i^ c , m ) J m +i(rep' a sin0) 
x cos(m0)} (2.101) 

where we have use the relation, given by 

cos(o; ± (3) = cos a cos (3 =f sin a sin /3 

and 

sin(o; ± (3) = sin a cos (3 ± cos a sin (3 
Substituting Eq.(2.58)~(2.61) into the Eq.(2.101)and (2.102), we have 

oo 

Afii = —j4:7ia 2 fi Go(R) cos9 V] i?on^ 2ri+ 5(Ksin6>)(/tsin#)~2 + 27ra 2 /ioGo(-R) cos ^ 

m=0 

oo oo 

x £ J m_1 £ {Ann[j m+2 n-i(«sin0) + J m+2n+ | (k sin 0)] (« sin 0)"5 

m=l n=0 

+5 mn [j m+2n+ i(/tsin0) - J m+2n+ |(resin0)](resin0)-i}cos(m0) 
+27ra 2 G'o( J R) cos# 



CXJ CXJ 

X E -T^ E {Cmn[j m+ 2n-\(KS™6) + J m+2n+ | (« sin 6)] (k sin 0)"3 
m=l n=0 

+Ann[4 i+ 2n+i( Ksin #) ~ J m+2n+§ ( K sin d )] ( K sin ~ * } Sm(m0) 

oo 

= -Ka 2 fi Go(R) E Con</ 2n+ | (/tsin^)(/tsin6 l )~2 + 2ira 2 G (R) 



(2.102) 



n=0 

■oo 



XE^ 'Ej^mn^+an-I^Sin^- J m+2n+|(« sin ^)](« sin ^) * 
m=l n=0 

+£ mn [j m+2n+ i(resin0) + J m+2n+ |(/tsin0)](/tsin0)-i}sin(m0) 
+27ra 2 G'o(i?) 



x 



oo oo 
E r 1 E {Cmn 



n=0 



^ m +2n-i( K sin0) - J m+2n+ s(K sin 9)](k sin ^ ± 



+D mn [j m+2n+ i(K sin 6) + J m+2n+ | (resin 0)] (resin 6) ^}cos(m0) (2.103) 
where we have again used the closure relation for the Hankel transform and the property 



of the delta function. 
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Next we evaluate A d z and given in Eq.(2.34) and (2.35), directly by applying the 
stationary phase method of integration. These integral can be written in the form, given 
by 

i nt = / o °° p(OJm(paO ex P [-Ve - K^Jr 1 ^ 

where we assume that -P(£) is slowly varying function. To perform this integration asymp- 
totically, we use the integral representation for the Bessel function given by 

jtn „+7r 

J m (paO = 7T exp(-j£p a cos a - jma)da 

Now transforming the cylindrical coordinate variable (p a , z a ) into the polar coordinate 
variables (R a ,9) through p a = R a sin9 and z a = R a cos9. And the integration variable £ 
can be change into by £ = /tsin/3. Then the integral I nt changes into 

-^nt = 7; — / exp(— jnRa sin/3 sin 9 cos a — jma)da 

/" ~ cos 3 

x / P(Ksml3)exp[-jKR a cosl3cos9}—^dl3 (2.104) 
jc sin p 

where the contour C is running along (0,0) — > (§,0) — >■ (§,00) in the complex /3-plane. 
Stationary points are located at (ao,po), which satisfied the equations 

sin 6> cos cos a — cos & sin Ad — 0, nR a sin # sin /3 sin a — m = 

when the value of nR a is sufficiently large, a may be set to 0, so that the approximate 
stationary points are given by 

«o = 0, 8o = 9 

Application of the standard process of the method yields the result, 

Ira = exph— — tt) P(K$\n9)—^- (2.105) 

v 2 ' nR a sm t/ 

Now applying this formula to the vector potential given in Eq.(2.34) and (2.35), we have 

A d g (r) = jpoa 2 -^ J -—J2B nJ 2n+ i^sm9)(Ksm9) * 

R sin 9 n=0 2 
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2 exp(-jK,R) cos 6 ~ rr . 

2.K sm& m=1 ra=0 

0) + D mn r~ n (Ksin0)] sin(m0)} 



+ 



r exp(-j/t J R) ~ . . 

p 2^ Con^n+f (« Sin 0)(« Sill 0) 

+eoa 2 Z eXp( ~f^ £ f;{[C mn S-„(«sin 

2jK m=l n=0 L L 

-A mn S~ n (/«sin6') + £ mn r+ n («;sm0) sinm^} 



cos(m0) 
(2.106) 



+ 



cosmcf) 

(2.107) 



In the far field only the and components of the E and H fields are dominant. Although 
the radial components are not necessarily zero, they are negligible compared to the and 
<j) components. Thus for far field observation, we have 



Eg = —juAg = ju sin 9A Z 

H e = -jcoFe = ju sin 6F Z = -YqE 6 

As = Z sin 9F Z 



(2.108) 
(2.109) 
(2.110) 



or 



where 
De(9,<j>) 



Eg 
Ea, 



exp(-jR) 
a D e {p,(p) 

exp(-jR) 



R 



(2.111) 
(2.112) 



—Z na cos 6 B 0n J 2n+ 5(n sin0)(/tsin0) 2 

n=0 

g oo oo 

+jy/ca cos J2 J m+1 J2{ [A mn ^ n (K sin 6) - B mn T^ n (K sin 6) 



m=l 



n=0 



+ [C„ m H+ n (/tsin0) + D mn r mn (nsin6 

oo ^ 

Z na C 0n J 2n+ 3(K sin 9) (k, sin 9)~ ? 



sin mc 



cos mcf) 
(2.113) 



n=0 
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17 OO OO 



oo oo 



1+ 

van 



(/tsin^) cosmcf) 



m=l n=0 



+ -A mn S mn (Ksin6>) + -B mn r+ n («sin0) sinm^} 



(2.114) 



These are the expression for the far-field pattern diffracted by a perfectly conducting cir- 
cular disk. Also, if we use EQ.(2.109)~(2.111) we find that Eq.(2.107) and (2.108) agree 
completely with Eq.(2.106). Therefore we will used Eq.(2.106) directly to find the far-field 
pattern diffracted by a circular hole in perfectly conducting plate. 

2.9 The Expressions for the Fields Diffracted by a Cir- 
cular Hole in a perfectly Conducting Plate 

The diffraction of electromagnetic plane wave by a circular hole in a perfectly conducting 
plane is a complementary problem of the scattering by a disk and the solution is obtained 
directly by using the result of the disk problem via Babinets's principle. 

2.10 Electric Filed Distribution 

From Eq.(2.38)~(2.41), the expression for the electromagnetic filed can be express in 
matrix form, given by 





-1 



(2.115) 
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H pc,m(Pa) 

H/f)s,m (Pa) 
Hps,m (Pa) 

H<j>c,m{Pa) 



/•OD r 

= / \H + (ZPa) 

Jo L 

r°° r 

= ±Y / \H-(Sp a 
Jo 1 



= ±Yn 



(2.116) 
(2.117) 
(2.118) 



where the kernel matrices [if ± (^p a )] are defined by 

'j'JfiPa) ±g/m(fr») 
±§- a Jm{tPa) JL(tPa) 



Using the boundary condition, the dual integral in this case can be written as 

pc,m(0 

-Ey>s,m(0 
Eps^mifC) 

-E</»c,m(0. 



/ [H-{iPa) 

Jo 1 



Z<% = Pa >l 
Z<% = Pa >l 



(2.119) 
(2.120) 



Y 



Y, 



f 

Jo 

f 

Jo 



H-{iPa) 



jVF^7?flfcm(0(«e) _1 







= 


_H<f>s,m(Pa) _ 






Hps,m(.Pa) 


= 


_Hlc,m\Pa) _ 





where H^ cm and -£/^, m represent the cosm0 and sinm0 parts of the incident wave H l p , 
respectively, and the same is true for H l ^ c m and . The expression for these factors are 
given by 



H; s ( P a) 
HIM 

HUpa) 



-jYqEx cos o e m j m J' m (K,p a sin 0„) 

jT E 2 e m j m ^-J m (np a sin 6> ) 

«p a sin 6*o 

jY E 2 e m j m J' m (K Pa sin 0„) 



jYb-Ei cos # e m y 



«p a sin 9 



sin^o) 



(2.123) 
(2.124) 
(2.125) 
(2.126) 
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The aperture electric field can be expanded in a manner similar to the disk problem. It 
is noted that (E p , E$) satisfy the vector Helmholtz equation V 2 E + /t 2 E = in circular 
cylindrical coordinates. Furthermore, (Ep^E^) have the property E p ~ (1 — p 2 a )~^ and 
E^ ~ (1 — p 2 )^ near the edge of the hole. By taking into these facts and using the property 
of the Weber-Schafheitlin's discontinuous integral, we set 



Epc t m(Pa) 



E, 



n=0 

oo 

pn.m(Pa) \CmnF mn (Pa) ~\~ D *mnG rrmiP a) 

n=0 

oo 

(Pa) = Y,[- A ™nF~ n (Pa) + B mn G+ n (Pa) 



n=0 

oo 



Efa m {Pa) — \pmnF mn (p a ) + D mn Gl^ n (pa) 



n=0 



(2.127) 
(2.128) 
(2.129) 
(2.130) 



where 



FL(Pa) 

G mn(Pa) 



J\m-l\(VPa)J\ m -l\+2n+±(v) ± ■Wl(*7Pa)4rH-2n+§ 7 7j 7 7 5<fo 7 ( 2 ' 131 ) 

/ [J\m-l\{VPa)J\ m -l\+2n+l(V) ± Jm+l(VPa)J m +2n+%V\V * d V (2-132) 
J 

where A mn ~ _D mn are the expansion coefficients and are to be determined from the remain- 
ing boundary condition that the tangential components of the magnetic field are continuous 
on the aperture. The corresponding spectrum function can be derived by applying the vec- 
tor Hankel transform, given by 

F pc,m{,Pa) 



Padpa 
Padp a 



(2.133) 
(2.134) 



F pc,m(£) 

_-E^s,m(0J [E^s^Pa) _ 

E ps ,m(0 (°°r ."I F pS:tn (p a ) 

Ea. (?) Jo 1 J frt 1 

Substituting Eq.(2.128)~(2.131) into Eq.(2.134) and (2.135) and perform the integration, 
then the spectrum function for the current density can be determined. The result is given 
by 



F pc ,m(£,) ^ [^■mn" mn (0 F mn T mn (£) 
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n=0 



(2.135) 



oo 

n=0 

oo 

Ep S ,m{Q = ^2[Cmn^ mn (0 + ^mnTmniO 
n=0 

oo 



n=0 



(2.136) 
(2.137) 
(2.138) 



2 E A 0n J 2n+ s(OC h E<fs,o(0 = "2 E ^onJ 2n+ |(Or § (2.139) 

n=0 n=0 

oo oo 

2 E Con-Wflr* ^c,o(0 = -2 E D, n J 2n Mi)C l (2-140) 



for m > 1 and 

Epc,o(Q ~~ 
E ps ,o(0 = 

n=0 " n=0 

for m — 0. From Eq.(2.116) and (2.117) we can express the spectral functions / cm (0 ~ 
<7sm(£) of the vector potentials in terms of those of the aperture distribution. This mean 
that the surface magnetic field can be expressed in term of the spectrum functions of the 
surface electric field. 



2.11 Derivation of the Expansion Coefficients 



From Eq.(2.116) and (2.117), the spectrum function / cm (£) ~ flw(0 ma Y t> e expressed 
in term of E(£), that is 



/cm(0£ j\f^ ^ Ep C ,m{£) 9sm(^)£, -^>s,m(0 



K 



1 



(2.141) 
(2.142) 



Substituting the above relations into Eq.(2.122) and (2.123), the continuities of the 
tangential components of the magnetic field are written as follows 

jy/( 2 -K 2 



Jo L 



3 ^Z^^PS,m\0 



H pc,m\Pa) 



= Pa <l 



(2.143) 
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f r 

Y \H + {i Pa ) 



Eq.(2.144) and (2.145) are projected into funct 



= Pa < 1 (2.144) 



ion 



. space with element v™{p 2 a ) for H p and 
u™(pa) for H^, then the matrix equation for the expansion coefficients A mn ~ D mn . The 



result are given by 



K' 



^m+2n-§(0 + ^m+2n+§(£) 



<^ +2P+ |(0 - k + 3 ) ■w w -|(o]de 



oo (if 



1 r 00 V£ 2 - /« 2 



—Km / 
Jo 



1 



/■oo 1 r I r I 1 

7 ^ ~ fgl _ ^ ^m+2n-±(0 — ^ m+ 2n+§(0 ^m+2p+§(0 + ^m+2p+§(0 

(0 - « + 3) (0 de 

^m+2n+§(0 + ^m+2ra+|(0 



'm+2n+|(0 ^m+2ri+§(£) 



-| 3 

3 ) J m +2p+l ( K sin #o)J sin 6» )~ 5 



= j2Ei cos # j m [a™ J m+2p+ 1 (k sin O ) - (a™ + 3) J m+2p+ | (k sin O )] (k sin 0„) 

m = l,2,3,..., p = 0,l,2,3,... 



(2.145) 



r 777 roc I r -, r , 

j E Cmn "J / V^-^ J m + 2„-l(0 + ^ m+ 2n + |(0] [^ + 2p-i(0 + ^ + 2p + |(0]^ 
n =0 ^ " 



/•oo ]_ r 

^^/q ^2 _ K 2 [ J m+2n-\{i) + ^m+2n+§(0 

oo 



°° { Tfl f°° ^/£ 2 — r 

-j E ^mn < — / -2 [J m+ 2n+\ (0 - ^ m +2r t +f (0 

„=0 I ^ J ° s 

a m J 

u, J, 



(0-K m + i)J m+2p+f (0 ^ 



^m+2p+i(0 + J m+2p+|(0j^ 



-ft 



/■oo J r 

'' _ K 2 [ J m+2n+\ (0 + ^m+2n+§ (0 

= — j2_Ei cos 9 j m m J m+2p _i (ftsin# ) + 



l+2p+ |(/tsin0 o ^ 
m = 1, 2, 3, • • • , p = 



_ i 

2 



ft sin #o 
',1,2,3,..., 



(2.146) 



/<=() 



1 /-oo - 



ft JO 



^m+2n-i(0 + ^m+2n+§ 
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f°° 1 r n 



°° (if 00 ^2 - ft 2 r ] 

+J X! -^mn S - / J m+2n+ l(£) ~ ^ m +2n+|(0 

n=0 I K J0 S 



(0 



—Km 



f°° i r 

_/ ^ 2 ^ /p _ [^m+2n+| (0 + ^ m +2n+| (0 



= -j2E 2 j m m[j m+2p+ i{Ksm6 ) + J m+2p+ |(/t sin# 

m = 1, 2, 3, • • • 



3 

P = 0,l,2,3,..., 

(2.147) 



00 (m f°° I r 1 r 1 

- J A mn \— \j ¥ ~ ^ J m +2n- \ (0 + ^ m +2n+f (0 J m +2p-\ (0 + ^ m +2p+f (0 ^ 

,/ (I I K JO - - 

"« /f 7|^[ J m+2n-§(0 " ^ + 2n + §(0] [<^ + 2 P -|(0 " « + l)<Wh§ (0] 

°° f 771 Z" 00 ^/^ 2 — ft 2 r -1 r -1 

+J E S ™M ~ / T J m+2n+\ (0 - J m+2n+§ (0 J m+ 2p-\ (0 + ^ m +2p+| (0 ^ 



/■°° 1 r 1 T m m 1 ] 

JO £\/ £ 2 — ft ? ^ m+2n+ ^^ ^"i+2n+f (0 ° ; ™^m+2p-|(0 ~~ I 01 ™ + l)^m+2p+§(0 ^| 

= -j2£ 2 j m [a™J m+2p _i(/tsin# ) - « + l)J m+ ^s(Ksin0o)](«sin0o)~5 

m = l,2,3, p = 0, 1,2,3, 



£A)n" / 



1 r ^2 



K 



J2U + I(0\-PJ2 P+ 1(0 + (P+ 1.5)J 2p+| (Ok£ 



r 13 

= -Ei cos 6>o -pJ 2p+ i (resin^o) + (p + 1.5) J 2p+ s («sin O ) (/tsin# )~ 5 



n=0 



o V£ 2 - /« 2 



J 2n+ a(0 -pJ 2p _i(0 + (P + 0.5)J 2p+ a(0 d£ 



2 P +; 
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— pJ 2p _i (ftsin# ) + (p + 0.5) J 2p+ 3(/tsin^o) (ftsin^o) 



P = 0,l,2,3, 



(2.149) 



2.12 Far Field Expression 



Far filed is obtained by applying the formula given in Eq.(2.106), to the expression of 
the vector potential, the result is given by 

2 V exp(-jKR) 1 



A*(p a ,(f>,z a ) = j2j2 a 2 Y - 



A 0n J 2n+ 3 (k sin 9) (k sin 9) * 



+fi a Y 
+ 



2 V <>Xp(- //.•/.') 



n=0 

oo oo 



sin# 

oo 

E ^ E {[A mn S- n ( K sine) - B mn r+ n (Ksm9) 

1 



m=l n=0 

D mn T^ n (K sin 9) 



smm 



sin^ 



cos m0 
(2.150) 



, . 9 exp(— ikR) cos9 



n=0 

oo oo 



+e a 



R sin 9 
2 exp(-j/ti?) 

m=l n=0 



R 



+ 



Ef" 1 E {[C mn H+ B («8ine) + L> mn r- n (/,sm^) 

cos 



-A mn S+„(/t sin 0) + B mn T mn (K sin 0) 



smm 



sin# 



cos mcf) 
(2.151) 



In the far region we have again the relation given by 

Eg = —juAg = jui sin 9A Z 
H e = -juF e = -Y E 4) 
Aa, = Z sin 9F Z 



(2.152) 
(2.153) 
(2.154) 



These are the expression for the far-field pattern diffracted by a circular hole in a perfectly 
conducting plate, when the observation point lie in the far region. 
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Chapter 3 



Numerical Computation 



The numerical computation has been carried out for the far-filed pattern diffracted by 
a perfectly conducting disk for normal incidence. The results are then compared with the 
physical optics (PO) approximation solution. The main problem in computational work 
is the numerical solution of the infinite integrals to obtained the expansion coefficient. 
Here we have transform these integrals into infinite series which are more convenient for 
numerical computation. Moreover, the numerical results for these integrals can be also 
obtained directly by using the numerical integration via the Gauss-Legendre quadrature, 
which is a simple check for the validity of these series solution. 



3.1 Physical Optics Approximation Solution 

We derive the approximate solution by using the physical optics to compare it with 
the results computed by the present method. We consider here the scattering by disk 
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illuminated by a plane wave. The induced current is given by 



J = 2n x H l = 2Y E 2 i x + E 1 cos6 i y exp(jK,xsm9 ) 



(3.1) 



The field produced by this current in the far region is derived from Eq.(2.87) as follows. 



2YoE 2 p Go(R) / exp jKp a sin 6 cos 0' + jnp a sin 9 cos(0 — 0') 
J s 

jKp' a Q cos(0' - $) p' a dp' a d(f>' 



dS 



= 2Y E 2 a 2 p G (R) / / exp 
Jo Jo 



AnY e 2 a 2 p G (R) C J (Kp' a Q) p' a dp' a d<f>' 
Jo 



where we set 



sin #o + sm cos = Ocos<3>, sin 6 sin = 9 sin <E> 

= yj (sin 6 + sin 6 cos 0) 2 + (sin 6 1 sin 0) 2 

The above integral can be readily carried out and we have 

2Ji(k6) 



A, = 27rFo^ 2 a>oGo(i?) : 



/t0 



Similarly we have 



A y = 2%Y E 1 a 2 cos9 p Go(R) 



2J 1 {kQ) 



Hence from Eq.(2.109), we have 



E e ~ — jcoA g = —2-kclk cos 6G (R) (e 2 cos + E 1 cos # s m 0) 

E<t> ~ ~~ J w ^0 — —j27raK,G {R) (—E 2 sin + E x cos # cos 0) — ~7T~~ 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



(3.6) 
(3.7) 



These are the expression for the electric field, which is obtained by using the physical optics 
method and can be used for the computational purpose, to compare the result with those 
obtained by the Kobayashi potential method. 
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3.2 Series Solution of the Integral G(a,(3]K) 



Consider an evaluation of the integral defined by 

G( a ,p; K ) = /f ^=fk (3-8) 

This integral was first evaluated by Nomura and Katsura. Here we have derive the series 
solution of this integral in a different way. That is, we define Eq.(3.8) as the limiting value 
of the integral given by 

G(a,#K) = lim r Ja{ ^^ -eM-\le - « 2 *K (3.9) 

In the first step, we drive the series representation for larger value of v(v > 1), then the 
expressions converted into contour integral and derive the result which is valid for smaller 
value of v(v < 1). Using the integral representation for the product of the Bessel function 
given by 

MOMO = 2 I' W2£cos0)cos([Ai-i/]0)d0 (3.10) 

7T JO 

And other integral formula for the Bessel functions given by 

j a+ M^e) = — r( „A/ f , u (Z™0) a+p+2t <it K<°) (3.ii) 

Zttj J-joo-€ 1 (a + p + t + 1J 

Eq.(3.8) is transformed into 

1 r+joo-t r(—t)dt r% 
G(a,/3-K) = — r ,„\i t , u / cos(^ +2t cos[(«-/3)^^ 

n z ] J-joo-f L(a + p + t + l) Jo 



x 



roo I / 

J f^'expl-^ - K*v]d£ (3.12) 



The integral with respect to 9 and £ may be carried out with the results 

7T 

K x = 1 2 (cos 6) a+ ? +2t cos[(a - 0)9\dB 
Jo 

nT(a + p + 2t + l) 



2 a +P+ 2t + 1 T(a + t + l)T(/3 + t + l) 
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(3.13) 



K 



7T. 



r 



-(a + /3)+t+ o 



2 [-Yi (a+0)+t (Kv) - jJi (a+ p)(Kv)] (3.14) 



( f )- 2+t 



Eq.(3.14) is valid for t> > 1. Substituting these results into Eq.(3.12) we have 

1 r +jco+i r(-t) T(a + f3 + 2t + l)T[l(a + /3) +t+ \] 



G(a, f3; k) 



A^pij J-joo-t T(a + (3 + t + l) 



T(a + t + l)T(/3 + t + 1) 



( ™L] 2+t 



(3.15) 



The Neumann function in the above equation is replace by the relation used as its definition, 
which is given by 

1 



Yi 



(a+f3)+t 



(kv) 



sin{[|(a + /3) +t]ir} 



x 



cos{(-(o! + 0) + t)ir}ji {a+ p )+t (Kv) - J_ { i {a+p)+t] (Kv)] (3.16) 



Using the above relation, G(a, (3; k) can be split into three parts, given by 

1 r+i°°-Z r/a + /3 



h = 

h = 

h = 

Q(t) = 



1 



Q(t) CSC ( 



J 



nv)dt 



Ay/ll J-joo-£ 
1 /-+.j'oo-£ 



Q(t)Ji 



(a+/3)+t 



(nv)dt 



^V^j J-joo-t 



Q(t) cot[(^y^ +t)rr 



Ji 



(a+j3)+t 



(nv)dt 



r(-t)r(a + f3 + 2t + l)r[|(a + £) + i + \] 



2±( ^_\ \{a+P)+t 

T(a + f3 + t + l)r(a + t + l)r(/3 + t + 1) K 2v J 



(3.17) 
(3.18) 
(3.19) 
(3.20) 



In the limit \t\ — > oo, the integral of l\ approaches to exp[2t \n(2/v )], while those of I2 
and J 3 approach to exp[— tlnn]. Therefore, we may close the contour of I\ in the left half 
plane for v < 2, and those of I 2 and J 3 in the right half plane. Since the evaluation of the 
integral depends on the indices a and j3, we consider the following cases, 
(a) Evaluation of the Integral I\ 

The integrand of the integral I\ has the singularities 
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(a) simple pole at t = -t (i = 1, 2, 3, • • • , ^) 

(b) double poles at t = _^±|±1 _ £ {£ = o, 1, 2, 3, • • • , ) 

(c) double poles at t = - Q+/3 + 1+£ (£ = 1, 3, 5, • • • , ) 

As it may be shown readily that the contributions from the double poles given in (b) and 
(c) vanish as v — > 0. Therefore, the contribution from only the simple poles give the result 
for h will be 

- §(a+0) r(£)r(a + ^ _ 2 £ + l)r(|(a + 



1 2 £ ^ r(a + /3-£+l)r(a-£+l)r(/3-£+l) 

Considering that — t is integer, we have 

te(y) J -^ )+ ,M = r( i (Q + p)- e + i) (3 - 22) 

Changing the index of summation from £ to p = — £, we finally obtain the expression 
of the series h 

i {a+ ^ 1)/2 r(±( a + 0)-p)r(2 P + i)r(p + ±) 



2A P =o r(|(a + /3)+p + l)r(|(a-/3)+p+l) 

x (3 23) 



(b) Evaluation of the Integral J 2 

The integrand of the integral I 2 has simple poles at t — p(p = 0, 1, 2, , • • • , ) and the result 
becomes the same form as the previous case. 



/ -v^y- r(q + /3 + 2p + l)r[|(q + /3)+p+|] 

2 " 3 ^h* T (P + X ) r ( a + P + P+ IS + /5 + 1) 



1 //cxa+/3+2p 

X r(/3 + +p + l)r[|(a + /3) + p + |] ^2 i ( } 
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(c) Evaluation of the Integral J 3 



This integral has simple at t — p + |(p = 0, 1, 2, • • • , ) and the result is given by 

- » (-i)pr(o! + p + 2 P + 2)r[|(o! + p)+p + i] 



I = -— 

3 2 ^ o r(a + (3 + p+l)T(a + p+l)r(P + p+l) 



x 



1 Q^ P+i (3 25) 



r(p + |)r[|(a + /3 + l)+p+l]V2- 



3.3 Series Solution of the Integral G 2 (a,/3;ft) 



Considering an evaluation of the integral defined by 



G 2 (a,fi,K) = J q eVF ^dZ (3.26) 
we defined Eq.(3.26) as the limiting value of the integral, given by 

G 2 (a,#K) = Urn jf 00 exp[-^ - ^t;]dg (3.27) 
Using Eq.(3.10) and (3.11), we have the above equation can be written as 

G 2 (a,fl«) = r°°"' w T^TTn / f cos"+^cos[(tt-/Wfl 

7r^j J-joo-£ L {a + p + t + I) Jo 

x ^_ g ^ +M " 2 exph^ 2 -^]^ (3.28) 

the integral with respect to and £ may be carried out i.e Eq.(3.13) and (3.14), the above 
equation becomes 

i r+j°o-s r(—t) 



T(a + f3 + t + l) 
T(a + f3 + 2t + l)r[|(a + 0) + * - §] («/2) a+ ^+ 2 *- 2 



x 



r(a + t + l)T(/3 + t + l) (/cu/2)( a +fl/2+t-i 
x -Fi^+t^H - jii( a+ffl+t -i(CT)]* (3.29) 
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The Neumann function in the above equation is replaced by using the relation, given in 
Eq.(3.16), G2(a, (3; k) can be splitted into three parts, given by 

1 r+joo-Z r /Cfc 4~ (3 \ 

1 r+joo-£ 

h = 7= / Q(t)Ji (a+w+( .iH* ( 3 - 31 ) 

1 f+joo-^ r /Q , _|_ , -i 

h = ^( t ) COt i(^^+ t - 1 ) 7r J J K^)+*-lM rft ( 3 - 32 ) 

nm r(-t)r(q + g + 2t + i)r[|(a + g) + 1 - \] (K , h(a+0)+t ^ ( . 

w; - r(a + /3 + t + l)r(a + t + l)r(/3 + t + l) W [ } 

In the limit \t\ — >■ oo, the integral of ii approaches to exp[2t ln(2/t>), while those of I2 and 
I 3 approach exp[— t £ v]. Therefore, we may close the contour of Ii in the left half plane 
for v < 2, and those of I 2 and 7 3 in the right half plane. Since the evaluation of the integral 
depends on the indices a and (3, we consider the following cases, 
(a) Evaluation of the Integral Ii 

The integrand of the integral ii has the singularities 

(a) simple poles at * = - \ -p (p = 0, 1, 2, 3, • • • , c ^ ± §= 1 ) 

(b) double poles at t = - q± §± 1 (p = 0, 1, 2, • • • , ) 

(c) double poles at t = - Q+/3 + 1+p (p = 1, 3, 5, • • • , ) 

As in the case (a), the contributions from the double poles vanish as v — > 0. The result is 

1 (a+ ^ 3)/2 r[|(a + p)-p- i]r(2p + 3)r(p + §) 



8 A p ^ T[\(a + P)+ V + 2]T[\(a-P)+ V + 2] 
T\\((3-a)+p + 2]T(p+l)(~2) (3 ' 34) 



x 



(b) Evaluation of the Integral J 2 

The integrand of the integral I2 has simple poles at t — p{p = 0, 1, 2, • • • , ) and the result 

becomes the same form as the previous one, that is 

■iyr(a + /3 + 2p + i)r[|(« + p) +p - \] 



3 8 ^ T{p+l)T{a + P + p + l)r{a+p+l) 
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1 /K\ a+/3+2p~2 

T{(3 + p + l)T[l{a + (3)+p] 



(3.35) 



(c) Evaluation of the Integral I 3 

The integrand of this integral has simple pole at t — p + \{p — 0, 1, 2, • • • , ) and the result 
is given by 

(-iyr(a + P + 2p + 2)T[l(a + P)+p] 



E 



x 



8 ^{a + (5 + p+l)Y{a + p+\)Y{(5 + p+l) 



3.4 Series Solution of the Integral K(a, (3] k) 

K(a, (3; k) is defined by 

K(a,frK) = lim jT y/? - K 2 J ^)M0 exp[-7e - « 2 *K (3.37) 
This integral may also be carried out as in G(a, (3; k). Here transforming K(a, (3; k) into 

H-i°°-« r(-*)r(« + p + 2t + 1) 



r 

K(a,p;K) = limlTrj y 



ioo-{ r(a + /5 + t + l)r(a + t + l)r(/3 + t + l)2°+/ ? + 2 *+ 1 

/•oo , 

X 



/ o °° V^ 2 - /c 2 r +/3+2i - 2 exp(-^ 2 - * 2 v)dt (3.38) 
Using the relation 

f°° 1 ( /rr«\ a+/3+2t //{n 2 /£n a+/3+2t+2l 



= 4 io TF^l^ -(2) (2) jexp^-^K 

^ r +/3+2t | T[§(a + /3) +t+ |] r x 



(Ku/2)3( a+ ^ +t 

r[l( a + i 9) + t-|] f , 

( ^ /2) i (a+ /3) +f -i {- y ic«+«+*-i(««) - jWiW) 
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K(a, (3; k) may be split into three parts 



K(a,p;K) = h + I 2 + I 3 

1 r+joo-Z 

h = -rr^-.l P(t){Q A (t)-Q A (t-l)}dt 
1 r+j°°-£ 

h = -— P(t){Q B (t)-Q B (t-l)}dt 

7T J-7'oo-f 



4^ J-joo-€ 

h = -j^ r°°^ P(t){Qc(t) - Qc(t - l)}dt 



(3.40) 



where 



r(-t)r(g + /3 + 2t + l) //s:xa+/3 + 2t 
U r(a + /3 + t+ l)r(a + t+ l)r(/3 + t+ 1) W l ' ' 

Q A(t) = r(^ +i+ i)(^)- , " +e+2 ' > c S c(^±4±^,)j_, (t , + , +21) ( K „)(3.42) 

0,(0 = r(^ + t + i)(^)- ( ° +W2 "j_, ( „ + , +2t) M (3.43) 

Qc(0 = r (^+ t+ 2)(y) cot(-(a + /3 + 2t)7r) 

X^(a+)8+2t)(«v) ( 3 - 44 ) 

The method of solving the above integrals are the same as discussed in the previous sections. 
Here we have show the result, that is 

= i r[i(a + /?)]r(2p+i)r(p-§) 

1 " V^F p ^ r(p+l)r[i(a + /3)+p + l]r[|(a + /3)+p + l] 

1 /K\2p 

X 



(=)* (3.45) 



r[|(a + /3) + p+l] v 2 

^ - (-l)pr(a + /3 + 2p + l)r[§(a + P)+p-\\ 
3 A 2.^ 



X 



4 p ^ r(p + l)r(a + /3 + p + l)r[|(a + (3) + p + 1] 

K\a+/3+2p 



(3.46) 



r(a + p + l)r(/3+p+l)V2 

~ (-iyT(a + /3 + 2p + 2)T[l(a + /3)+p] 
4 ^ r(a + p + p + f)r(p + §)r[|(a + 0) + p + §] 

1 /K\ a+/3+2p+l 



X 



r(a + p+f)r(/3 + p+f)V2 
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It is noted that K(a, (3; k) is related to G(a, (3; k) and G 2 (ck, (3; k) by the relation given by 

K(oi,P;k) = G(a,f3;K) -G 2 (a,f3;K) (3.48) 



3.5 Numerical Integration 



Numerical integration is the approximate computation of an integral using numerical 
techniques. The numerical computation of an integral is sometime called quadrature. The 
most straightforward numerical integration technique uses the Newton-Cotes formulas, 
which approximate a function tabulated at a sequence of regularly spaced intervals by 
various degree polynomials. If the endpoints are tabulated, then the 2- and 3-point for- 
mulas are called the trapezoidal rule and the Simpson's rule, respectively. The 5-point 
formula is called Boole's rule. A generalization of the trapezoidal rule is Romberg integra- 
tion, which can yield accurate results for many fewer function evaluations. 

If the function are know analytically instead of being tabulated at equally spaced in- 
tervals, the best numerical method of integration is call Gaussian quadrature. By picking 
the abscissas at which to evaluate the function, Gaussian quadrature produces the post 
accurate approximation possible. However, given the speed of modern computers, the ad- 
ditional complication of Gaussian quadrature formalism after makes it less desirable that 
simply brute-force calculating twice as many points on a regular grid. Here we have com- 
puted the integrals by using the Gauss-Legendre quadrature, which is simply a check for 
the series solution of the integrals. 
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3.5.1 Numerical Integration of G(a,(3;K,) 



The integral G(a, (3; k) can be also computed numerically by using the Gauss-Legendre 
quadrature. Here we set 

G(a,/3;K) = G re (a,(3; k) - jG im (a,(3; k) (3.49) 

The imaginary part of G(a,(3; k) is computed from 

G im {a^-K) = j K J_ J a (OJff(0^ 

JO y K — 

= 2 w L= J a { K - V 2 )Mk - V 2 )d V (3.50) 

The real part is split into three parts as 

G re (a,/3;«) = f°° l^ 2 - FUOMOdt = G 1 + G 2 + G 3 (3.51) 

where 

Gi = / - T ==J a (0M0dZ = 2 - 7 ===J a (K + r ] 2 )MK + r ] 2 )dr ] (3.52) 
Jk V K — q JO v + 7/ 

G2 = £t^W j " K)j ^ kr (3 ' 53) 

r°° 1 z* 00 1 / k 2 \ 

g 3 = / o -j=p=^j a (t)Mt)dz = ] x _(i + _)j a (o^(Ode 

1/1 1 [K 2 / , ,\-n (a - /3)7r 1 ai - 6i . (a - /3)7r 

" + 3^ iy - (° 2 + h ~ aibl )\ I cos —r~ + n^2x*- sm —r~ 

11 . ( n v a + (3 + 1 x 01 + 61 1 / a + p + 1 



where 



sm(2X tt) - cos(2X X) (3.54) 



4a 2 -1 (4a 2 -l)(4a 2 -9) 
01 = — 8" ' ° 2 = 128 
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3.5.2 Numerical Integration of the Integral G2(a,(3;K,) 



We set again 

G 2 (a,(3;n) = G 2re (a, (3; k) - jG 2im (a, (3; k) (3.56) 

The imaginary part of G is computed from 

r K 1 

G 2 im(a, p; k) = / c2 J a {£)Jp{xi)<K, 

JO \J K — 

= 1 r J a (K - V 2 )Mk - v 2 )dr, (3.57) 

The real part is split into three parts as 

[00 I 

G 2re (a,P;K) = / J a (Z)MO<% = g 2i + G 22 + G 23 (3.58) 

where 



2k I 

= U i^wv^ JjK+l ' 1)MK + '' 1)dl ' (3 ' 59) 

If 1 1 r/t 2 / 7 \ 1 1 (a-p)n lai-fei . (a - f})ir 

- ^4X^ Sm ( 2X 2 )-^^X^ CO < 2X 2—") (3 - 61) 

where ai,6i and a 2 ,6 2 are defined in Eq.(3.55). 



3.5.3 Numerical Integration of the Integral K(a, (3; k) 
We set 

K(a,P;K) = K re (a,(3;n) + jK im (a, j3; k) (3.62) 
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The imaginary part of K(a, 8; k) is computed from 



/oo . / l~2 _ C2 



v 2 ) 

The real part is split into three parts as 



2 i 2K ~ 2 t v 2 Ja(K ~ V 2 )M* ~ V 2 )dv (3.63) 



K re (a,/3; K ) = / V , 2 q J«(fl M<x)d£ = + ^2 + ^3 (3.64) 



where 



^^a(O^R 



£ 2 

lo (k + t] 2 

rX ^2 _ K 2 



2/ , " I V 2 UK + V 2 )MK + V 2 )dv (3.65) 



J3k c; 
/■oo — K 2 

k 3 = j x V ^ 2 uoum 

roo 1 , if2 



I r 1 1 r/t 2 / , , m -1 (a — 8)tv 1 ai — 61 fa — /JW 

II a + /3 + lx 1 ai + 61 / v a + /3 + l 
sin ' ' " ' 



vr 2X2 ""^ 2 V ~ ^ C ° S F 2 V (3 ' 67) 

To verify the validity of the series expression of G(a,(3;K), G 2 (a, (3;k) and K(a,f3;n) 
which are derived in the previous section, we perform numerical computation and the result 
are compared with those obtained by direct numerical integration. The agreement is fairly 
good. 
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Chapter 4 



Numerical Discussion 



With the formulation developed in the previous section, we have preformed the computa- 
tion of the diffracted field for a disk with na = 5.0, kcl = 9.0, kcl = 12.0 and kcl = 16.0. For 
all cases computed, the thickness of the perfectly conducting disk is considered as infinites- 
imal. In order to check the effectiveness of the present method based on the Kobayashi 
potential, the numerical data for all the cases mention above are compared with the results 
obtained using the physical optics (PO) method. 



4.1 Computations of Matrix Elements and Radiation 
Pattern 

A first step in obtaining numerical results of the physical quantities is to compute the 
matrix elements defined in Eq.(2. ) and (2. ). These are infinite integral and can be 
computed by using the series solution of the infinite integral. From this, we obtained 
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the expansion coefficient by using the matrix inversion. Once the numerical result for the 
expansion coefficients are obtained, the radiation patterns are computed from Eq(2.80) and 
(2.84). The numerical results of the far-filed pattern diffracted by a perfectly conducting 
plate for normal incidence are shown in figures below. The ordinate a denotes the power 
pattern or the differential scattering cross section. To verify the validity of the present 
computation, we show the result produces by the physical optics (PO) in these figures. 

4.2 Conclusion 

We have formulated the electromagnetic field diffracted by circular conducting disk and 
circular hole in the conducting plate, when plane wave impinges on the obstacle, by using 
the method of the Kobayashi potential. We have write the expression for the vector po- 
tential by using the Helmholtz equation, in the form of the spectrum function. From these 
vector potential the corresponding electromagnetic field are derived. These expression are 
then written in the matrix form for the sack of simplicity. Finally the field expressions 
are expanded in the form of double series and each summand satisfies desired edge condi- 
tions as well as a part of other required boundary conditions. The mathematical formulation 
involves dual integral equations derived from the potential integral and the boundary con- 
dition on the plane where the disk or hole is located. The dual integral has been solved 
by using the Weber-Schafheitlin's discontinuous integral and the the orthogonal property 
of the Jacobi's polynomials. From this we obtained the weighting function in form of the 
matrix equations. Expansion coefficients are determined from the solution of the matrix 
equation and the simple series expansions for the matrix elements are derived. The matrix 
equation involves the infinite integral, which has been transform into infinite series expan- 
sion. These series expansion are more convenient for numerical computation. 

We derive the far field expressions by two different ways. One method is to evaluate 
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the field radiated from the current density induced on the disk and the second method 
is to evaluate the expression of the vector potential directly by applying the stationary 
phase method of integration. We presented the numerical results of the far diffracted field 
pattern by normal incidence, which is fairly agree with physical optics method (PO). The 
main problem in computational work is the numerical solution of the infinite integral which 
involve in the matrix equation for the expansion coefficient. Here we present the infinite 
series solution of the these infinite integral. Furthermore the the direct numerical solution 
of these integral are also presented for the compression with the series solution. The present 
method promises applicability to a wide class of problem such as, circular hole in the thick 
conducting plate, flanged circular resonator, and so on. 
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Figure 4.1: Far- field pattern diffracted by a perfectly conducting plate for normal incidence 
(0 O = 0°, = 0°). Disk is k = «a = 5. 
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Figure 4.2: Far-field pattern diffracted by a perfectly conducting plate for normal incidence 
(0 O = 0°, = 0°). Disk is k = na = 9. 
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Figure 4.3: Far-field pattern diffracted by a perfectly conducting plate for normal incidence 
(0 O = 0°, <p = 0°). Disk is k = na = 12. 
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Figure 4.4: Far-field pattern diffracted by a perfectly conducting plate for normal incidence 
(0 O = 0°, <p = 0°). Disk is k = na = 16. 
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Appendix A 

Vector Hankel Transform 



The vector Hankel transform is the generalization of the conventional Hankel transform 
and it transform a vector function from one space to a vector function in another space. 
This was introduced by Chew and Kong [41], [42] to analyse the disk resonator and antenna. 
Here we will reproduce the proof by assuming the transform pair 



h{pa 

fl{pa 

HO 



f 

Jo 



/ [H^p-a) 
Jo 1 



m). 

fl(Pa) 
f2(Pa) 



Padpa 



where 



J'miiPa) 
J21 



Jm{iPa) 



±^J m {iPa) J'miiPa) 



If Eq.(Al) is true, it implies 

fl(Pa) 
f2(Pa) 



(Al) 



(A2) 



(A3) 



(A4) 
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Eq. ( A4) can be proved by interchanging the order of integration. The diagonal element of 
ff±(fro)] [ff^fr'j] become 

777,^ 777, d 

£J'm{£Pa)J' m {£Pa)+t J J m (fro) Jm (fro) = (4+1 (fro) 4+1 (fro) + 7" 37 4 (fro) ■/ to (frl) 

sPaPa PaPa a ? 

Using the closure relation of the Hankel transform, given by 

<5(Po - Pi) 



/■oo 

/ fr m +l(fr a )J m +l(fr' a )^ 



Po 



the first and second term evaluates to zero. Similarly, the off diagonal elements can be 
written as 

Tfl III III 

-JLttPaVrnWa) + -■Ufra)4»(&'a) = r J ™(fro)4(fr a ) 

Pa Pa PaPa 

which integrate to give zero. Consequently Eq.(A4) becomes 



Pa^Po 

which complete the proof. 



£df/ P' a dp'a ff ± (fra) ^(fr'J 
JO L J L 



Pa 

Pa 



/ 2 (P'J 





7i(po)" 




/l(Po) 




./ 2 (p a ). 







(A5) 
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Appendix A 

Weber-Schafheitlin's Discontinuous 
Integral 



The Weber-Schafheitlin's discontinuous integral is defined by 

W{fi,u,X;r)= dfi (A.l) 

JO u 



To solve the above integral, we used the relations 



' 1 ~\z/+2m 



jf x ) = y (-i) m (a.2) 



and 



roc 

/ exp[-at]J„(&t)*' t ~ 1 dt 
Jo 



^ 1 m\T u + m + l) Jo i [ J 



m=0 

oo 



m\T(v + m+ 1) 



T( ir (l^ +2m r(/x + i/ + 2m) 

; m!r(z/ + m + l) a"+ I/ + 2m 1 ' ; 



Hence we can write 
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W(n, u,\;r) 



r"T[±(fi, + v - A + l)] 

2 A r[|(- /U + z/ + A + i)]r(z/ + i) 

r/i + Z/-A + l pL — U — X + l 

xF 



< r < 1 



m(/i + z,-A + i)] 



2V- A + 1 r[|(ji - 1/ + A + l)]T(u + 1) 
/i + v — A + 1 —fj, + i/ — A + 1 



1 1 



r > 1 



(A.4) 



2 7 2 

If there exists a relation /i — v + A = —2m — l(m = 0, 1, 2, • • • , ) among the parameters, 
the function becomes W(/i, is, A; r) = for r > 1 because of the properties of the gamma 
function|r(— n)\ —¥ oo (n: positive integer). For r in the range < r < 1, W(/j, : v, A; r) 
takes finite value. 
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Appendix A 



Expansion of the Bessel Function by 
Jacobi's Polynomials 



A.l Projection of J m ±i(€y/x) into v™(x) 

We project J' m {i\/x) and j^J m (^y/x) into functional space with element v™{x). To 
evaluate the components of each element , we need the following relation 



h = 



j\^ m+1 \l-x)h m ^(^K(x)dx 
" V8(2p + m + j)T(q + m+ j) Vh# 1>ro 

r(9 + i)r(m) £§ 
v^r(m + i)r(p+|) ^p+m+i V8r(m + i)r(p + |) ^ P+m +f 

(2p + m+ §)r(p + m) £§ (2p + m + |)r(p + m + 1) £§ 

v^r(m + i)r( P + |) r / 2 2„+ m +i(0 / 3x^ 2p + m +§(0 

( 2, +m+ i)r ( p +m+ i) — a — 



(A.l) 
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V8(2p + m + |)r(g + m + f) ^g+m+f (0 
^ T(g + l)r(m + 2) 



jm+l,m 

q,p 



8p7(m + l)r(m + |) ^ 2p+m +i(0 



tf 3 



(2p + m + |)r(p + m+ 1) 

m + |)r(m + l)r(m + l)T(p + f) ^ 2p+m +f (0 
(2p + m + |)r(p + m + 1) ; J 



/8(p 



Jo 



v^r(m + i)r(p + |j 

(2p + m+ |)r(p + m+l) 



~P -3 + ( P + 

£ 2 



3\ ^2p+m+ 

2 J 



v^r(m + i)r(m + f, 

772 

(2p + m+ |)r(p + m+ 1) 



2p+m+i(0 ^2p+m+§(0 
£ 2 ^ 2 



v^r(m + i)r(p + |) 

(2p + m+ §)r(p + m + l) 



x 



• t in t 2 ) L \P t "i -r i ; 

(2p + m) ^li® _ (2p + m + 3 )%±i® ' 

£ 2 £ 2 J 



have used the relations 



jm—l,m 



jm—l,m 

n,n 



f x m {l-x)^^ l \x)^{x)dx 

J 

r(m)r(m + l)r(n + 2)r(n + §)r(2n + m + |) 
F(2n + m + |)r(n + m + l)r(m + n + §) 

; x m (l-x)^™" 1 (a;)z/™(a;)dx 

T(m)T(m + l)r(n + l)r(n + §)r(2n + m + \) 
r(2n + m + f)r(n + m)r(m + n + I) 



, m+ \)T{n + m)T{m + n + \) 

JO 

r(m + l)r(m + 2n)r(n + l)V(n + |)r(m + 2n + §) 



™+l,m _ /" 
n.n / 
JO 



T(2n + m + |)r(n + m + l)T(m + 
m+1 {l - x)^v™ +1 {x)v™{x)dx 



JO 

T(m + l)r(m + 2)r(n + l)r(n + §)r(2n + m + |) 
r(2n + m + |)r(n + m + l)T(m + n + |) 
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Projection of J m ±i(£y/x) into u™(x) 



oo 



^ r( g + i)r(m) ^ 
v^r(m + i)r(p+§) v^r(m + i)r(p + §; 

(2p + m + + m) 



^2p+m+§ 



(0 



y/2T( 



m + i)r(p + |; 



(2p 
fl 



|)r(p + m + l) 



(p + m) 



(2p + m + |)r(p + m + 1) £: 

■^2p+m-i(w / l\^2p+m+|(0 



£ 5 



(A.9) 



.|(m+l) ( 



''(1 - X)~5 J m+1 (£y/x)u™(x)dx 



V2(2p + m+l)T(q + m + l) ^ g+m+ |(0 lm 

h r(g + i)r(m + 2) e § 

v^ P r(m + i)r(p + 1) ^p+m-iCO v^fo + m + |)r(m + i)r(p + 1) ^ 2p+m +§(0 

(2p + m + |)r(p + m + 1) ^1 (2p + m+ |)r(p + m + 1) 



v^ P r(m + i)r(p + ^) 

(2p + m + |)r(p + m + l) 



-p- 



+ (P + 



m + 



2p+m+f 



(0 



(A.10) 



j 1 x^ m+1 \l - 
Jo 



m 



x , 



J m (ev^)<(a:)d 



r(m + i)r(p 



>/2(2p + 77l+ §)r(p + m+l) 



rn 



J. 



2p+m- 



X 



;(0 ^ ^2p+m+| 



(0 



r(m + i)r(p+|) 



"3 



v/2(2p + m+ |)r(p + m+ l)2m(m + 2p+\) 
1 _ Ji m (^) u ™{x)dx 



(0 



(A.H) 



r( m + i)r(p + |) 



^(2^ + 777 + |)r(p + 777 + 1) 



(2p + 777 J j-^ (2p + 777 



1)- 



J, 



2p+m+ 



(A.12) 



K 



m—l,m 
n+l,n 



j 1 x m (l 
Jo 



X) 



x)dx 



r(777)r(m + i)r(n + 2)r(n + j)r(2w + m+J) 

r(2n + 777 + |)r(n + m + l)r(777 + 77 + |) 



(A.13) 
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J x m (l-x)- 1 *u™- 1 (x)u™(x)dx 

T{m)T{m + l)T{n + l)T{n + §)r(2ra + m - \) 
T(2n + m + |)T(n + m)T(m + n — |) 

JO 

r(m + l)r(m + 2)r(n + l)r(n + ^)r(2n + m-\) 
r(2n + m + |)r(n + m + l)r(m + n+\) 

f 1 x m+1 {\ - x)-tu™ +1 (x)u™(x)dx 
Jo 

T(m + l)r(m + 2)r(n + l)r(n + |)r(m + 2n + |) 
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